Fidelity for multimode thermal squeezed states 



Gh.-S. Paraoanu* 

Department of Physics, University of Illinois at Urbana- Champaign, 1110 W. Green St., Urbana, IL 61801, USA; e-mail: 

paraoanu@physics.uiuc.edu 

Horia Scutaru 

Department of Theoretical Physics, Institute of Atomic Physics, PO BOX MG-6, R-76900 Bucharest-Magurele, Romania; 

e-mail: scutaru@theorl . theory. nipne. ro 



In the theory of quantum transmission of information the 
concept of fidehty plays a fundamental role. An important 
class of channels, which can be experimentally realized in 
quantum optics, is that of Gaussian quantum channels. In 
this work we present a general formula for fidelity in the case 
of two arbitrary Gaussian states. From this formula one can 
get a previous result [Q], for the case of a single mode; or, one 
can apply it to obtain a closed compact expression for mul- 
timode thermal states. The concept of fidelity used in this 
paper is the standard one pl-W. It can be defined by 



max |('!/'i|'(/'2)| 

\i>i>,\ip2> 



where >, i = 1, 2 are purifications of the density matrices 
Pi- 

03.65.Bz, 42.50.Dv, 89.70.+C 



I. INTRODUCTION 

Within recent years, the quantum theory of informa- 
tion, an extension of the classical theory of information 
to the quantum realm, has emerged as a fascinating re- 
search field. A great deal of effort has been devoted to 
the issue of transmitting a state through noisy quantum 
channels despite the quantum-mechanical uncertainties 
in our knowledge about that state. This is a different 
problem from the classical situation, where the states are 
mutually exclusive and the input system may remain in 
its initial state; in the quantum case the states are density 
operators on a Hilbert space and the non-cloning theo- 
rem precludes the input system, in general, to retain 
its original state. Moreover, for a noisy quantum chan- 
nel the state is subject to a decoherence process, due to 
the interaction with an external environment, which fur- 
ther decreases the reliability of information processing. 
Thus, a fundamental problem is to extend the classical 
encoding and decoding procedures to quantum channels 
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and to define an upper limit (the channel capacity) to 
the amount of quantum information that can be trans- 
mitted with an arbitrary high fidelity. Recently, work |^ 
has been focused on quantum channels which use Gaus- 
sian (or quasi- free, (tJ^) states for the transmission of 
information. In this first issue to be raised is how 

to calculate the fidelity given two mixed Gaussian states; 
partial answers, depending on the particular type of mix- 
tures under consideration have already been given in the 
literature §,1,0. In this paper, we will give a general 
formula for the fidelity of two quasi-free states, and show 
that the previous results can be obtained as particular 
cases. The formula can also be applied to some interest- 
ing cases, such as two-mode systems 

Let pi and p2 two density operators which describe two 
mixed states. The transition probability P{pi, P2) has to 
satisfy the following natural axioms: 

1. P{pi,p2) < 1 and P(pi,p2) = 1 if and only if pi — 

2. P(pi,p2) = P(p2,Pi); 

3. If p\ is a pure state, p\ = >< ■0i| then 
P{pi,P2) =< ipi\p2\il'i >; 

4. P{pi,p2) is invariant under unitary transforma- 
tions on the state space; 

5. P {pi\a, P2\a) > P{pi,P2) for any complete subal- 
gebra of observables A; 

6. P(pi ® (Ti,p2 (E) (J2) = P{pi,P2)P{(Jl,(72)- 

Uhlmann's transition probability for mixed states |3| 

P{pi,P2)=[Tr{^P2^f']' (1) 

satisfies properties 1-6. The fidelity is defined by 
F{pi,p2) = P{pi,p2). A detailed analysis of the struc- 
ture of the transition probability was hampered by the 
factors containing square roots in Eq. (^. Due to tech- 
nical difficulties in the computation of fidelities, few con- 
crete examples of analytic calculations are known. Un- 
til recently, all the results were obtained only for finite- 
dimensional Hilbert spaces [|2|-Q . The first results in an 
infinite-dimensional Hilbert space were recently obtained 
by Twamley M for the fidelity of two thermal squeezed 
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states and by Paraoanu and Scutaru [|0| for the case of 
two displaced thermal states. In Scutaru has devel- 
oped another calculation method which allowed getting 
the result for the case of two displaced thermal squeezed 
states in a coordinate-independent form. 

Let {E, a) be a phase space i.e. a vector space with a 
symplcctic structure cr. Then the commutation relations 
on (E, a) acting in a Hilbert space Ti. are defined by a 
continuous family of unitary operators {V{u),u £ E} on 
H which satisfy the Weyl relations 

V{u)V{v) ^exp^^^a{u,v)^V{u + v). (2) 

Hence the family {V{tu), — oo<i<cx3}fora fixed u E E 
is a group of unitary operators. 
Then by the Stone theorem 



V{u) — exp{ii?(M)}, 



(3) 



where R{u) is a selfadjoint operator. From the Weyl 
relations we have 

exp{«ii?(M)} exp{isi?(u)} — 

exp{its(T(u, v)} exp{isi?(w)} exp{iii?(u)}. 



By differentiation and taking t — s — one obtains 
[R{u),R{v)] = ~ia{u,v)I. 



(4) 



The operators {R{u),u G E} are called cannonical ob- 
servables. 

The phase space E is of even real dimension 
2n and there exist in E symplectic bases of vec- 
tors {sj, fj}j^i^,,,^n, i-e. reference systems such that 
(j{ej,ek) = (j{fj, fk) = and (j{ej, fk) ^ -(j{fk,ej) = 
5jk, j,k = The coordinates rj^) of a vector 

M e -E in a symplectic basis, u = -I- rj^ fj), 

are called symplectic coordinates. The well known coor- 
dinate and momentum operators are defined by Qk = 
R{fk) and Pk = -R(efc) for k = l,2,...,n. Then the 
canonical observables R{u) are linear combinations of 
the above defined coordinate and momentum operators: 

There is a one-to-one correspondence between the sym- 
plectic bases and the linear operators J on E defined 
by Jgfc = -fk and J fk = Ck, k = The es- 

sential properties of these operators are: (t( Jtt, u) > 0, 
a{Ju,v) + aiu^Jv) = {u,v G E and = — /, where 
/ denotes the identity operator on E). Such operators 
are called complex structures. In the following we shall 
use the matriceal notations with u € E as column vec- 
tors. Then (t{u,v) = Jv and the scalar product is 
given by a{Ju,v) = u^v,u,v G E. A linear operator S 
on E is called a symplectic operator if S'^JS — J. When 
S" is a symplectic operator then 5*-^ and are also 
symplectic operators. The group of all symplectic oper- 
ators Sp{E,a) is called the symplectic group of {E,a). 



The Lie algebra of Sp{E, a) is denoted by sp{E, a) and 
its elements are operators R on E with the property: 
[JR]^ — JR- Hence an operator R on E belongs to 
Sp{E, a) n spiE, a) iff i?^ = _ j. if j and K are two 
complex structures, there exists a symplectic transfor- 
mation S such that J = S^^KS. For any symplectic 
operator S we can define a new system of Weyl opera- 
tors {V{Su); u e E}. Then from a well known result on 
the unicity of the the systems of Weyl operator up to a 
unitary equivalence it follows that there exists a unitary 
operator U{S) on H such that V[Su) = U{S)'^V{u)U{S). 

For any nuclear operator O onH. one defines the char- 
acteristic function 



CFu{0) =TrOV{u), ugE. 



(5) 



We give the properties of the characteristic function 
which are important in the following 0: 

1. CFo(O) = TrO; 

2. CFu [V{v)*OV{v)] = CFu [Oexpicr(i;,u)]; 

3. CFuiOi02) = 

4. CFsuiO) = CFu{U{S)OU{S)^). 

II. MULTIMODE THERMAL SQUEEZED STATES 

The multimode thermal squeezed states are defined by 
the density operators p whose characteristic functions are 
Gaussians pP,|i 



(6) 



CFu{p) = exp "I ~^u^Au 



where Ais a2nx2n positive definite matrix, called corre- 
lation matrix. From the last property of the characterisic 
function, enumerated above, it follows that: 



A 



U{S)pU{S)^ 



S AqS 



(7) 



Because the correlation matrix A is positive definite it 
follows IqJlql that there exists S € Sp{E, cr) such that 



A^S'VS 



(8) 



where I? = ^ ^ ^ J and 13 > / is a diagonal n x n 

matrix. The most general real symplectic transformation 
S e Sp{E, cr) has pUT^] the following structure: 



where 



M 



S^OMO 



M 




(9) 



(10) 
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and O, O are symplectic and orthogonal {O^O = /) 
operators, and where M is a diagonal nx n matrix. Var- 
ious particular kinds of such matrices are obtained taking 
O, O , V or M to he equal or proportional to the cor- 
responding identity operator. A pure squeezed state is 
obtained when T) = I. If this condition is not satisfied, 
the state is a mixed state called thermal squeezed state 
10. When A4 = I there is no squeezing and the corre- 
spondig states are pure coherent states or thermal coher- 
ent states. All these states have correlations between the 
different modes produced by the orthogonal symplectic 
operators O and O . As a consequence the most general 
form of a correlation matrix A is given by: 



A 



o"^MO'^VOMO' 



(11) 



From the property 3 of the characteristic function we 
have for two density operators pi an p2 



CF„(pip2) 
1 

exp 



det 



A1+A2 



[A2 - {A2 - iJ)iA, + A2)-\A2 + iJ)] 



When pi = p2 we have 



CFu{p ) = [det A) 2 exp ^-^u 



1 T f A- JA-^J 



(12) 



A state p is pure iff = p. Then from the equality 
CFuip^) = CFu{p) it follows that a Gaussian state is 
pure iff 



A = -JA-^J, 



(13) 



i.e. a Gaussian state is pure iff J A £ Sp{E,a). Anal- 
ogously, for a mixed state < p. Then CFu[p^) < 

CFu{p) and as a consequence "^—^ — - > A. Hence for 
any Gaussian state the correlation matrix A must satisfy 
the foUoving restriction 



A < - JA-\J. 



(14) 



III. THE CHARACTERISTIC FUNCTION OF 
THE SQUARE ROOT OF A DENSITY MATRIX 

Let us suppose that the characteristic function of the 
Hilbert-Schmidt operator ^/p of a Gaussian state is, up 
to a numerical factor, also a Gaussian function with the 
correlation matrix <&(A), 



CF.,{^) = if exp |-iu^<i>(A)u| . 
Then from the equality p = ^/py^ we obtain 



(15) 



(det($(A)))- 



exp 



exp 



$(A) - J$(A)-iJ 



'-u'^Au 
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Hence 



and 



= v/rfei$(A), 



$(yl) - J^A)-^J = 2A. 
The last equation has the solution 

$(A) = A{I + v// + (JA)-2). 



(16) 



(17) 



(18) 



(19) 



This is a new proof of a result obtained in 0|. The ad- 
vantage of this proof is given by the fact that it does 
not require the choice of a special basis in i?. If we take 

a symplectic basis in E such that ^ ^ g ^ then 

JV = VJ and from this equation and from the equa- 
tions A = S^VS, JS^ = S-i J it follows that (JA)-^ = 
S-^V-^S. Hence / + (JA)-^ = S-\l - V-^)S and 



AyJl+{JA)-^ = S^{VV^-I)S 



(20) 



IV. THE GENERAL FORMULA FOR THE 
FIDELITY OF GAUSSIAN STATES 

The fidelity F(pi, P2) for two density operators pi and 
P2 is defined by 



F{Pl,P2) = Tr dy/p{p2y/Pl 



(21) 



As we have seen in section II the characteristic function 
of a product of operators whose characteristic functions 
are Gaussians is also a Gaussian. In section HI we have 
obtained a simple formula for the characteristic function 
of the square root of a density operator whose character- 
istic function is a Gaussian. Hence we can find a simple 
formula for the characteristic function of the operator 
JTiP-iJpi- 



(22) 



CF,{^^P2^i) = Vlexp|-iz^O^ 



where 



= det^{Ai)-'det \^ 

, .'A2 + <P{Ai)~U 
det 



/$(Ai) + A2 



(23) 
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where U = {A2 - iJ)($(^i) + ^2)"^(^2 + iJ), and 
O = - ($(Ai) - iJ)[A2 + $(Ai) - 

Then applying the result of the preceeding section we can 
obtain the characteristic function of \J yfp\P2\fp\-, 



CF, 



P1P2VP1 



[Ldet^{0)Y exp { --z' $(0)^ 



(24) 



From this formula and the property 1 of the characteristic 
function we obtain 



We remark that 



def^iO) = detOdet I + ./T+{JO)^ 



(25) 



(26) 



In order to simplify the formula for fidelity we observe 
that 



ti]k = TrpiPjPk = det 



A, + A, 



Aj + Ak- {Aj - iJ){A, + A^y^iAj + iJ) 



det 



and that ^123 = ^231 = ^312- If we take in this last identity 
<i>(Ai) instead of Ai then we obtain 



det 

det 



$(Ai) + ^2 



A2 + ^{Ai)-U 



det ( ^l±Al ) det<^{Ai). 



Hence we get 



L 



det 



Ai + ^2 



(27) 



It is not evident from this general formula that the 
properties 1-6 of the fidelity are valid. Let us consider 
the most simple one, namely the property F{p,p) — 1. 
In this case it is necessary to prove that $(0) = A. We 
can choose the complex structure J to commute with 
the correlation matrix A: J A = AJ. Then all opera- 
tions in the formula which gives O as a function of A 
and J can be performed and the result is: O — — 
and $(0) = A. The next property which we shall dis- 
cuss is the property 3 which in the case of Gaussian 
states becomes 0: F{pi,p2) = (^^4^)"^- We shall 



prove that O = Ai when pi is a pure state. First 
we remark that ^{Ai) = Ai and that there is a sym- 
plectic transformation such that Ai = S^S. Then 
O^S^{l- 2P_ [21 - 4P+A'-ip_] 2P+} S where 

= Id^, P_ = and X = {S'^)-^A2S-^+L Evi- 
dently P+ and P_ are an orthogonal decomposition of the 
unit operator: P^ = P+, = P_, P+P_ = P_P+ = 
and P-f -I- P_ — I. As a consequence of the orthogonality 
we obtain O = S'^S = A^. Then $(0) ^ $(Ai) = Ai = 
S^S and det{^{0)) = det{S'^)det{S) = 1. 



V. THE ONE MODE CASE 

In we have obtained an expression for the fidelity 
in the one mode case. This formula can be reobtained as 
a consequence of the above general formula. In the one 
mode case all matrices are 2x2 matrices. For a 2 x 2 
matrix O we have 



$(0) = eO, 



(28) 



where e = 1 + y/l - and det^{0) = {VdetO + 
Vdet O — 1)^. From these considerations it follows that 

2 



F{pi,P2) 



v/ctelpT+Aa) ( VdetO - \/detO - 1) 



(29) 



Thus it is sufficient to compute detO. We shall denote 
by V the product {detAi — l){detA2 — 1). After simple 
but long computations we obtain 



detO = 1 



det{Ai+A2) 
which gives the result of fll 

2 

F{Pl,P2) 



y/det{Ai+A2)+V-Vv' 



(30) 



(31) 



VI. MULTIMODE THERMAL STATES CASE 

In the case of two thermal states with correlation ma- 
trices Ai = Vi with i = 1, 2 we have AiJ = JAi, {i = 1, 2) 
and A1A2 = A2A1. Then $(A,) = A, + ^If^, 
(j = 1,2). Hence 



and 



$(0) 



O ^ {Ai + A2y\AiA2 + I) 



{A, + A2) 

iA,A2 + I)-ViA'i-I){Ai-I) 



(32) 



(33) 



Finally 



4 



F{pi,P2) = 



det 



\ \iA,A2 + I) - ^{A^ 
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which is the product of the fidelities of the corresponding 
one-mode thermal states Ml 



VII. CONCLUSIONS 



In this paper we have provided a general formula for 
the calculation of the fidelity of two Gaussian states. It 
is shown that, in the particular case of a single mode, 
this formula reproduces the results already known in the 
literature, and in the case of multimode thermal states it 
yields a compact expression with a direct interpretation. 
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